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COMPUTATION OF SCALAR FAR-FIELD PATTERNS OF 


LARGE-APERTURE ANTENNAS 
by Thomas A. O'Malley 
Lewis Research Center 

SUMMARY 

In computer programs used for evaluating the performance of high- gain antennas, 
efficient numerical methods for calculating the far-field patterns must be used since 
the majority of computer time and storage requirements may be attributed to this 
phase of the program. The numerical method most frequently used is the Fast Fourier 
Transform (FFT), which computes the far field as the Fourier transform of the field 
distribution in the antenna aperture. This report describes a new numerical method 
that in many applications is superior to the FFT in terms of reducing computer time 
and storage requirements. 

In the new method, the arbitrarily specified field in the aperture is represented by 
a Fourier series expansion. The far field is obtained by evaluating a double integral 
where the integrand contains the aperture function. When the Fourier series expres- 
sion is substituted into the integrand, the integral can be evaluated in closed form. 

The new method requires less computer storage than the FFT method. Which of 
the two methods requires less computer time depends on the number and location of 
far-field points. If the application requires far-field points only in the principal planes 
of the pattern, the new method is generally faster. If the application requires a full 
rectangular array of far-field points, the FFT is generally faster. A significant ad- 
vantage of the new method is the complete flexibility in choosing the number and loca- 
tion of far-field points. If the calculation of the field in the aperture is time consum- 
ing, the number of aperture grid points required becomes an important consideration. 
The number of aperture grid points required is approximately the same for the two 
methods. 



INTRODUCTION 


Reflector antennas and lens antennas whose dimensions are large compared to a 
wavelength are used extensively in applications requiring high gain. To determine the 
performance characteristics of these antennas, accurate calculation of the antenna 
pattern in the far field is required. One method of calculating the far- field pattern of 
such antennas is the aperture- field method (ref. 1). In this method, the field vectors 
E and H are specified in a planar aperture in the vicinity of the antenna. The far 
field can be obtained from vector diffraction equations involving integration over the 
aperture. These vector equations are complex and require considerable computational 
effort. In many antennas the field over the aperture is almost completely linearly po- 
larized. Only a small fraction of the energy is in the cross-polarization component of 
the field. By neglecting the cross-polarization components, the problem can be re- 
duced to a scalar diffraction problem, which requires considerably less computational 
effort. 

Scalar diffraction theory can be used to show that the far field is the two- 
dimensional Fourier transform of the field distribution in the aperture of the antenna 
(ref. 1). In computer programs used for evaluating antenna performance, efficient 
numerical methods for calculating the two-dimensional Fourier transform must be 
used since the majority of computer time and storage requirements may be attributed 
to this phase of the program. The most widely used method for calculating the far 
field is the Fast Fourier Transform (FFT), which, as its name implies, is consider- 
ably faster than previously used methods (ref. 2). 

This report presents a new numerical method for calculating from scalar diffrac- 
tion theory the far- field patterns of high- gain antennas. The method evaluates the far- 
field pattern as the two-dimensional Fourier transform of an arbitrarily specified field 
distribution in the antenna aperture. It is shown to be superior to the Fast Fourier 
Transform in many applications. 

The new numerical method is described in detail. A brief discussion of the Fast 
Fourier Transform follows. A comparison of the two numerical methods for three 
common aperture distributions is then given. Criteria for determining which of the 
two numerical methods is superior in a given application are developed. 


DESCRIPTION OF NEW NUMERICAL METHOD 

In this new numerical method, we assume that the antenna aperture lies in the x-y 
plane (fig. 1). In nearly all cases of interest, the aperture is either elliptical or cir- 
cular. Let R be the vector from the center of the aperture to a far- field point. Let 
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a be the angle between the z-axis and the projection of R into the x-z plane, and let 
/? be the angle between the z-axis and the projection of R into the y-z plane. Let the 
coordinates of a point in the aperture be x, y. The field over the aperture will be 
designated by F(x, y). A(x, y) will be the amplitude distribution, and <p(x, y) will be 
the phase distribution, so that 


F(x, y) = A(x, y)e^ (x > y) (1) 

The far field will be designated by E(q:, /3). From reference 3, the far field is given 
by 


E(o, 0) =ff F(x, y) e j< 27r/x >< x shl sin dx dy 


( 2 ) 


where the integration is over the aperture A and X is the wavelength. Let us define 

(3) 


k = — — sin a 
x \ 


k =— sin (3 
y X 


Equation (2) then becomes 


If 1 


^Vv.' F(x, y)e 
A 


^(k^-tkyy) 


dx dy 


Let us define the function u(x, y) as 

u(x, y) = F(x, y) inside A 
u(x, y) = 0 outside A 

Then equation (5) becomes 

-j(k x x+kyy) 


/» OO /» oo 
J—oo •/— oo 


Eflv k y ) = 


u(x, y)e 


dx dy 


(4) 


(5) 


( 6 ) 

(7) 


( 8 ) 


The far field expressed in this form is the two-dimensional Fourier transform of the 
aperture function. 
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Fourier Series Representation of Aperture Ill uminat ion 


In the numerical method developed herein, the aperture illumination function 
F(x, y) is represented by a Fourier series expansion. When the Fourier series expres- 
sion is substituted into equation (2), the integral for the far field can be evaluated in 
closed form. 

The aperture is assumed to be elliptical in shape. The elliptical aperture includes 
the circular aperture as a special case. Let us consider an extension of the aperture 
illumination function to a rectangular region enclosing the elliptical region, as in fig- 
ure 2. The enclosing rectangle will be called the Fourier rectangle. The semimajor 
and semiminor axes of the elliptical aperture are a and b, respectively. The semi- 
length and semiwidth of the Fourier rectangle are c and d, respectively. On this 
rectangle the aperture illumination can be approximated by a Fourier series S^x, y) 
of the form 


M M 

F(x,y)=»S M (x,y)s £ c rs e 

= -M s=-M 


J(r7rx/c) e j(STry/d) 


The Fourier coefficients c _ are given by 

rs 


c rs 


-A- f° / d F(x,y)e“j< r7rx/c >e-j< s,r y /d >dydx 

4cd 4 J-d 


(9) 


-M ^r^M; -M < s ^ M (10) 

We assume that the aperture illumination is specified on a set of grid points within the 
Fourier rectangle. We further assume that the grid points form a (2L + l)-by-(2L + l) 
array of equally spaced points, where L is some integer. The grid points (x., y^) are 
given by 

x. = -c +■ 2ci i = 0, . . ., 2L (11) 

1 2L + 1 


y =-d +— -M_ j = 0 , . . ., 2L (12) 

J 2L + 1 

To evaluate the Fourier coefficients in equation (10), we replace the double integral by 
the following double summation: 
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(13) 



-j(nrx./c) -j(S7iy./ d) 
F(x., yj )e e 1 


-M ^ r s M; -M s s s M 
To calculate the c__ in equation (13), one may, of course, evaluate the summations 

rb 

directly. Instead of doing this, we use a more efficient algorithm devised by Goetzel. 
For details of the algorithm, see reference 4. 

The values of c and d and the method of extending the aperture illumination 
function outside the elliptical aperture are chosen so as to obtain a good Fourier series 
approximation within the elliptical aperture. For grid points outside the elliptical 
aperture, a good approximation is not required because the far-field pattern is ob- 
tained by evaluating an integral over the aperture only. 

The Fourier series approximation S^-(x, y) has the periodicity property that 

s M (c, y) = S M ( _C ’ y) for a11 y 311(1 s m^ d > = s m* x ’ -d ) for a11 x - m g eneral > th e 

aperture illumination function does not have periodicity properties. Thus, if c = a 
and d = b, the Fourier series approximation may be poor in the vicinity of the points 
(a, 0), (-a, 0), (0, b), and (0, -b). As an example of this poor approximation, consider 
figure 3, which is a plot of the imaginary part of the aperture function and its Fourier 
series representation as a function of x, with y set equal to 0. The aperture function 
is 


F(x, y) = < 0. 3+0.7 



j(7r/2)(x/a)' ; 


(14) 


In figure 3, a = b= c = d=l, and M = 4. It is clear from figure 3 that the approxima- 
tion is poor at x = ±a. Figure 4 is the same plot but with c = 1. 1 a and d = 1. 1 b. 
The Fourier series approximation at x = ±a is significantly better. Furthermore, 
the approximation is significantly better throughout the entire interval. Numerical 
tests suggest that the ratio of c to a and the ratio of d to b should be between 1.05 
and 1. 10 to obtain good approximation throughout the entire aperture. 

The method of extending the aperture function outside the aperture should be done 
such that the function has the same periodicity properties as its Fourier series repre- 
sentation; that is, F(c, y) = F(-c, y) for all y and F(x, d) = F(x, -d) for all x. The 
method used is illustrated in figure 5. Let y be fixed at y = y^ and consider the row 
of grid points (x^, y.), i = 0, . . . , 2L + 1. For some value of i, let us say i = il, the 
grid point (x.p y.) lies just inside the aperture in the second quadrant. For some 
value of i, let us say i = i2, the grid point (x- 2 , y.) lies just inside the aperture in the 
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first quadrant. The extension of the function F(x, y) is such that 


F(x 0 , y.) = F(x 2L+1 , y.) = ± [f^, y.) + F(x- 2 , y.)] (15) 

The function F(x, y^) is extended by two linear segments: one extending from (Xq, y.) 
to (Xjp y.), and the other extending from (x i2 , yj) to O^l+P y j)» as s ^ own in figure 5. 
These extensions are performed for all such rows of grid points. The same procedure 
is then used to extend all columns of grid points. 


Evaluation of Diffraction Integral 


Substituting the Fourier series representation of F(x, y) into equation (2), we ob- 


tain 


/Yv v 

E(a,«= / / 2/ °rs e 

J J r=-M s=-M 


j(r7rx/c) e j(S7ry/a) e j(27r/X)(x sin ot+y sin p) dy ^ (16) 


In the appendix, the following integral is evaluated: 



e j(7x4 l5 y) dydx = 


2*abJ 1 V(ya) 2 + (6b)' 




(ya) 2 + (6b) 2 


Combining equations (16) and (17), we obtain 

M M 

E(o:, P) = 27rab 


I Z’" i tr 


r=-M S=-M 


where z"L is given by 
rs 


rs 


_ / r 7 ra f 27ra sin cA 2 + / S7rb 27rb sin /3 

\ c X / \ d X 


(17) 


(18) 


( 19 ) 


The computer time required to evaluate equation (18) is reduced by prestoring the 
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function f(x) = J 1 (v'x)/4/x for a set of values of x. After z^ g is calculated by equa- 
tion (19), the expression J-i (z„ 0 )/z„_ is evaluated by linear interpolation on the table. 

i rs rs 

Since the far field is expressed as a sum of Bessel functions, this numerical method 
will be called the Bessel Function method, or simply the BF method. 


FAST FOURIER TRANSFORM 

The Fast Fourier Transform (FFT), as its name implies, is a numerical method 
for calculating Fourier transforms that is computationally faster than previously used 
methods. Since the far field can be expressed as the two-dimensional Fourier trans- 
form of the aperture illumination, the FFT can be used to perform the computation. 
For details of the FFT algorithm, see reference 5. We discuss briefly here only 
those features of the FFT needed to make a comparison between the FFT method and 
the BF method. 

The FFT requires that the aperture illumination be specified on a rectangular 
array of N-by-N points. As in the discussion of the BF method, we call this rectangle 
the Fourier rectangle. From figure 6, the Fourier rectangle has a semilength of c 
and a semiwidth of d. Unlike the BF method, the FFT method requires that the aper- 
ture illumination be set to zero for those points that are outside the aperture. Also 
tinlike the BF method, the Fourier rectangle is large compared to the aperture. The 
reason for this is that the ratio of c to a (or d to b) determines the number of 
points per side lobe calculated by the FFT. Approximately c/ a (or d/b) points per 
side lobe will be calculated. 

The number of grid points lying within the aperture determines the number of 
lobes calculated by the FFT. if the number of grid points on the x-axis, extending 
from x=-a to x = a, is k, the FFT calculates approximately k lobes. Also the 
number of grid points lying within the aperture determines the accuracy of the far- 
field calculation. The accuracy increases as k increases. 

The Fourier- rectangle array of N-by-N points has the restriction that N must 
be a power of 2. For most applications, N will be either 64, 128, or 256. Another 
restriction is that the number of far-field (output) points must also form an N-by-N 

array, where N is the same as for the Fourier- rectangle (input) array. 

2 

The computer time for calculating the FFT is proportional to N log 9 N. Previ- 

3 " 

ously used methods had computer times proportional to N . For large values of N, 
the savings in computer time by using the FFT is substantial. 
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COMPARISON OF BF METHOD WITH F FT METHOD 


For either the FFT method or the BF method, the error at a far-field point is de- 
fined as the difference between the computed and true values at that point, divided by 
the peak value of the true pattern. By true value, we mean the exact value of the 
scalar diffraction integral of equation (2). The error criterion used in this report is 
that the error cannot exceed 0. 01 at any far-field point. This criterion will not yield 
good accuracy on a decibel scale for side lobes far removed from the main lobe. How- 
ever, far-out side lobes are not accurately calculated from the scalar diffraction inte- 
gral (ref. 6). Vector diffraction integration is required. The error criterion of 0. 01 
is acceptable for those portions of the far field that are accurately described by the 
scalar diffraction integral. 

For the FFT method, the accuracy of computing the scalar diffraction integral in- 
creases as the number of lobes included in the calculation increases. To satisfy the 
accuracy criterion of 0.01, approximately 16 lobes have to be included. And to obtain 
visual resolution of the side lobes, approximately 8 points per side lobe are needed. 
Thus, a 128- by- 128 array of points is required. 

For the BF method, the upper limit M in the far-field summation has to be 4 to 
satisfy the error criterion of 0. 01. For some aperture illuminations, a smaller value 
of M would suffice, as discussed later in this section. For calculating the Fourier 
coefficients c rg , a 17-by-17 array of points is sufficient for the Fourier rectangle. 

In the sample aperture illuminations to follow, c = 1. 05 a and d = 1. 05 b for the BF 
method. 

The first aperture illumination to be considered is given by 


A(x, y) = 0. 3 + 0. 7 



■\ 


> 


(p(x, y) = 0 


( 20 ) 


This illumination is typical for large- aperture antennas. We assume that the aperture 
is circular and thus a = b. For the BF method, M = 2 is sufficient to satisfy the 
error criterion. Figure 7 is a plot of the relative power lE(u, 0)/E(0, 0)1 in decibels 
as a function of a. The FFT method, the BF method with M = 2, and the BF method 
with M = 8 are shown. The BF method with M = 8 can be considered as the exact 
evaluation of the scalar diffraction integral since larger values of M do not produce 
significant changes. 

We now consider an aperture illumination having a quadratic, or square- law, 
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phase function. This phase function occurs in large- aperture antennas when the feed 
is defocused along the axis of symmetry. For the BF method, M = 4 is required to 
satisfy the error criterion. The aperture illumination is given by 


V' 


A(x, y) = 0. 3 + 0. 7 



<p(x, y) 



-\ 




( 21 ) 


Figure 8 shows the FFT method, the BF method with M = 4, and the BF method with 
M = 8. Again, the BF method with M = 8 can be considered as the exact evaluation 
of the scalar diffraction integral. 

Finally, we consider an aperture illumination having a cubic phase function. This 
phase function occurs in large- aperture antennas when the feed is moved off axis to 
tilt the beam. For the BF method, M = 4 is required to sa’tisfy the error criterion. 
The aperture illumination is given by 


A(x, y) = 0. 3 + 0. 7 



<P(x, y) 


7 r 
2 





( 22 ) 


Figure 9 shows the FFT method, the BF method with M = 4, and the BF method with 
M = 8. The latter can be considered as the exact evaluation of the scalar diffraction 
integral. 

In comparing the merits of the two methods, one should consider the computer 
storage and computer time requirements. The FFT requires an N-by-N array of 
complex numbers for input. Thus, 2N real numbers must be stored. The FFT uses 
the same storage space for calculating the N-by-N array of complex outputs. For 
N = 128, the storage requirement is approximately 32 000 words. The BF method re- 
quires 2000 words of storage for storing the table of the function J^(Vx)/x/x. For a 
K-by-K array of aperture points, an additional storage of 2K^ is required. For typ- 
ical values of K, this storage requirement will not exceed 1000 words. Finally, the 
output points must be stored. The maximum storage requirement would occur when a 
full rectangular array of N-by-N output points is desired. Since only the amplitude of 
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the far field is of interest, the storage requirement is N words. The total storage 

2 2 

requirement for the BF method is 2K + N + 2000 words. For any reasonable values 
of N and K, the BF method requires less storage than the 2N words required by 
the FFT. 

In discussing computer time requirements, there are two major parts of the com- 
putation to be considered. One part is the computation of the aperture illumination at 
the aperture grid points. For idealized aperture illuminations, such as those consid- 
ered in this report, this computer time is insignificant. But for illuminations calcu- 
lated from feed horn patterns, the computer time may be significant. For the two 
methods, the computer time is approximately proportional to the number of grid points 
within the aperture. In the examples used in this report, a 16-by-16 array of aperture 
points was required for the FFT method, and a 17- by- 17 array was required for the 
BF method. It may be concluded that the computer time for calculating aperture il- 
luminations is comparable for the two methods. 

The second major part of the computation is the calculation of the far field at the 

output points. For an N-by-N array of output points, the computer time requirement 

2 

for the FFT method is proportional to N loggN. Since loggN will be either 6, 7, 

or 8 in most applications, we may conclude that the computer time is approximately 

2 2 
proportional to N . For the BF method, the computer time is proportional to M J, 

where M is the upper limit in the double summation for the far field and J is the 

number of output points. 

To determine the computer time for calculating the far field at the output points, 

we use results obtained from the IBM 360/67 time- sharing system. These results may 

not be the same for another computer and, therefore, should be taken only as guide- 

2 

lines. For the FFT method, the computer time was 1. 2 N milliseconds. For the BF 
method, the computer time was 1. 6 M J milliseconds. 

In many cases, the far field needs to be known only in the two principal planes of 
the pattern. Suppose that N points per plane are required. Then the BF method re- 
quires 2N output points and the FFT method requires N output points. Of course, 
in the FFT method, we obtain more information than we need because a full rectangular 
array of output points is obtained. We compare the computer time needed for both 
methods for N of 64, 128, and 256, assuming M = 4 for the BF method, in the fol- 
lowing table: 
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Method 

Value of 

N 


64 

128 

256 


Computing time, sec 

FFT 

aX 

H 

21 

96 

BF 

3 

6 

13 


For such cases, it is clear that the BF method is faster than the FFT. 

For some applications, the full rectangular array of output points is required. 

2 

Then N output points are required for both methods. For this case the computer 
time for the BF method divided by the computer time for the FFT method is 1. 3 M . 
Clearly, the FFT method is faster for this case. 

A significant advantage of the BF method is the complete flexibility in choosing 
the number and location of output points. Recall that the FFT method is restricted to 
having the number of output points equal to the number of input points and that this 
number must be a power of 2. The BF method has no such restrictions. Also short- 
cuts can be taken with the BF method by taking advantage of symmetries in the pattern. 
For a circularly symmetric pattern, for example, output points are needed only for 
one-half of one principal plane of the pattern. 


SUMMARY OF RESULTS 

In computer programs used for evaluating the performance of high- gain antennas, 
efficient numerical methods for calculating the far- field patterns must be used since 
the majority of computer time and storage requirements may be attributed to this 
phase of the program. The numerical method most frequently used is the Fast Fourier 
Transform (FFT), which computes the far field as the Fourier transform of the field 
distribution in the antenna aperture. This report describes a new numerical method 
that in many applications is superior to the FFT in terms of reducing computer time 
and storage requirements. 

The new method requires less computer storage in all applications. Which of the 
two methods requires less computer time depends on the number and location of the 
far-field points. If the application requires far-field points only in the principal planes 
of the pattern, the new method is generally faster. If the application requires a full 
rectangular array of far-field points, the FFT method is generally faster. A signifi- 
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cant advantage of the new method is the complete flexibility in choosing the number and 
location of far-field points. If the calculation of the field in the aperture is time con- 
suming, the number of aperture grid points required becomes an important considera- 
tion. The number of aperture grid points required is approximately the same for the 
two methods. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 4, 1976, 

643-60. 
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APPENDIX - EVALUATION OF INTEGRAL 


The following integral is to be evaluated: 



(Al) 


(A2) 


The integration over y is readily done to obtain 



•/-a 


dx 


(A3) 


Since the imaginary part of the integrand of equation (A3) is an odd function of x, the 
imaginary part of I is zero. The trigonometric part of the integrand can be ex- 
pressed as 


sin 




(A4) 


Substituting equation (A4) into (A3), we obtain 



Making the substitution 


x = a sin 6 
dx = a cos 9 d 6 


(AS) 


(A6) 
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we obtain 


I - — 

35 


Z 7r/2 /* 71-/2 

^•(ya sin 0+6b cos 0) CQS e de _ f e j(ya sin 0-6b cos 0) cog Q de 

_ /2 J-iti 


■7r/2 


The following substitutions are made: 

ya sin 0 + 6b cos 0 = A sin(0 + 0 ^) 
ya sin 0 - 6b cos 0 = A sin(0 + a 2 ) 
where A, <r 1 , and cr 2 are given by 


= /(ya) 2 + (6b) 1 




The integral I now becomes 


1 = -^ 

36 


jAsin(0+t7,) r 

/ e cos 0 d0 - J e * cos 0 d0 

^-7r/2 •^■tt/2 


•7t/2 . 


jA sin(0-hr 2 ) 


(A7) 

(AS) 

<A9) 

(A10) 

(All) 

(A12) 


(A13) 


Since I is real and a/j6 is purely imaginary, the real parts of the two integrals in 
equation (A13) are zero. These two integrals, to be called 1^ and l^, are evaluated 
in the following equations: 

/’ 7r/2 jA sin(0-hj.) 

I. = / e 1 cos 0 d0 i = 1, 2 (A 14) 

A/2 
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,(7r/2)4tr. 

* 4 > 

- (7r/2) -NTj 


/•( ^ 

L = / e^ A sin e cos(e - tr.Jde i = 1, 2 

./-fir 


(A15) 


■ (7r/2)-+tj. 


I. = cos or. f e^ A sin e cos e de 

*^-(ir/2)+Oj 


+ sin cr. y 


(7t/2)+CT. 

(7r/2)-HT. 


e jA sin e sin £ dg j - lf 2 (A16) 


Since the real parts of Ij are zero, we obtain 
,(7r/2)-hr. 


= j cos a. / 

J-(-, 


sin(A sin e)cos e de 




.(7r/2)-Hr. 


r 

+ j sin a. / sin(A sin e)sin e de (A17) 

J- (it/ 2) +or. 

The integrands in equation (A17) are periodic, with period 7r, so we can change the 
limits of integration to range from -(l/2)7r to (1/2) ir. The first integrand in equation 
(A17) is odd, and the second integrand is even. We therefore have 


r */ 2 

I. = 2j sin u. J sin(A sin e)sin e de i = 1, 2 
From reference 7, we obtain 

I. = j7r sin or. J-^(A) i = 1, 2 
Returning now to the evaluation of I, we have 


a7rJ,(A) 

1 = . (I 1 " V = (sin °1 “ sin °2> 

j 6 6 


(A18) 


(A19) 


(A20) 
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From equations (All) and (A12), it follows that 



Using equations (A10), (A20), (A21), and (A22), we finally obtain 



(A21) 

(A22) 

(A23) 
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Figure 6. - Coordinate system for antenna aperture and Fourier rectangle - FFT 
method. 



Figure 7. - Far-field pattern for aperture distribution 1: A(x,y) = 0.3 + 0.7(1 - lx 2 +v 2 )/a 2 ]- 
<plx,y)*0. 


20 


Relative power, dB Relative power, dB 



Figure 8. - Far- field pattern for aperture distribution 2: Alx.y) s 0. 3 + 0.7 [1 - lx 2 +y 2 )/a 2 ] ; 
iplx.y) = (7t/2)[(x 2 + y 2 )/a 2 l . 



Figure 9. - Far-field pattern for aperture distribution 3: A(x,y) = 0. 3 + 0. 7 [ 1 - (x 2 + y 2 )/a 2 ]; 
<p(x,y) = (7r/2)(x/a) 3 . 


NASA- Langley, 1976 


E-8746 



NATIONAL AERONAUTICS AND SPACE ADM I N lSTRATlO^J 
WASHINGTON, D.C. 20546 


OFFICIAL BUSINESS 

PENALTY FOR PRIVATE USE *300 


SPECIALS FOURTH-CLASS RATE 
BOOK 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 
451 



596 001 Cl 0 D 760618 S00903DS 

DEPT OF the air force 

,.»o 

KIRTLAND AFB NM 87117 



POSTMASTER : 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute . . . to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



